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Abstract — Sensorless control of Permanent-Magnet Syn- 
chronous Motors (PMSM) at low velocity remains a challenging 
task. A now well-established method consists in injecting a high- 
frequency signal and use the rotor saliency, both geometric 
and magnetic-saturation induced. This paper proposes a clear 
and original analysis based on second-order averaging of how 
to recover the position information from signal injection; this 
analysis blends well with a general model of magnetic saturation. 
It also proposes a simple parametric model of the saturated 
PMSM, based on an energy function which simply encompasses 
saturation and cross-saturation effects. Experimental results on a 
surface-mounted PMSM and an interior magnet PMSM illustrate 
the relevance of the approach. 

Index Terms — Permanent-magnet synchronous motor, sensor- 
less position estimation, signal injection, magnetic saturation, 
energy-based modeling, averaging. 

I. Introduction 

PERMANENT-Magnet Synchronous Motors (PMSM) are 
widely used in industry. In the so-called "sensorless" 
mode of operation, the rotor position and velocity are not 
measured and the control law must make do with only 
current measurements. While sensorless control at medium 
to high velocities is well understood, with many reported 
control schemes and industrial products, sensorless control at 
low velocity remains a challenging task. The reason is that 
observability degenerates at zero velocity, causing a serious 
problem in the necessary rotor position estimation. 

A now well-established method to overcome this problem is 
to add some persistent excitation by injecting a high-frequency 
signal |1] and use the rotor saliency, whether geometric for 
Interior Permanent-Magnet machines or induced by main flux 
saturation for Surface Permanent-Magnet machines El- flOl . 
Signal injection is moreover considered as a standard building 
block in hybrid control schemes for complete drives operating 
from zero to full speed fTTll - |[T5ll . 

However to get a good position estimation under high- 
load condition it is important to take cross- saturation into 
account llT6ll-lf26ll . It is thus necessary to rely on a model 
of the saturated PMSM adapted to control purposes, i.e. rich 
enough to capture in particular cross-saturation but also simple 
enough to be used in real-time and to be easily identified in the 
field; see 1271 -1321 for references more or less in this spirit. 

The contribution of this paper, which builds on the prelimi- 
nary work [33], is twofold: on the one hand it proposes a clear 
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and original analysis based on second-order averaging of how 
to recover the position information from signal injection; this 
analysis can accommodate to any form of injected signals, 
e.g. square signals as in (34], and blends well with a general 
model of magnetic saturation including cross-saturation. On 
the other hand a simple parametric model of the saturated 
PMSM, well-adapted to control purposes, is introduced; it 
is based on an energy function which simply encompasses 
saturation and cross-saturation effects. 

The paper runs as follows: section [III presents the saturation 
model. In section [Till position estimation by signal injection 
is studied thanks to second-order averaging. Section [IVl is 
devoted to the estimation of the parameters entering the 
saturation model using once again signal injection and averag- 
ing. Finally section HV-Cl experimentally demonstrates on two 
kinds of motors (with interior magnets and surface-mounted 
magnets) the relevance of the approach and the necessity of 
considering saturation to correctly estimate the position. 

II. An energy-based model of the saturated PMSM 
A. Notations 

In the sequel we denote by xij := (xi,Xj) T the vector 
made from the real numbers X{ and xj, where ij can be dq, 
a/3 or j 5. We also define the matrices 



cos fi — sin fi 1 
v sin/i cos/i 

and we have the useful relation 



dM^ 
dfi 



JCM U 



and JC 







B. Energy-based model 

The model of a two-axis PMSM expressed in the syn- 
chronous d — q frame reads 



dt 
J duo 

ri 2 dt 
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with <pdq flux linkage due to the current; 



(A,0)^ 



(1) 
(2) 
(3) 
con- 



stant flux linkage due to the permanent magnet; Ud q impressed 
voltage and id q stator current; uj and 6 rotor (electrical) speed 
and position; R stator resistance; n number of pole pairs; J 
inertia moment and tl load torque. The physically impressed 
voltages are u a p := MgUdq while the physically measurable 
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currents are i a p := Mgidq- The current can be expressed 
in function of the flux linkage thanks to a suitable energy 
function H((j)d,(t>q) b Y 



i dq -± dq w dq) .- y d2U ^ d ^ q) 



(4) 



where d^H denotes the partial derivative w.r.t. the k th vari- 
able 03, E3; without loss of generality 7/(0,0) = 0. 
Such a relation between flux linkage and current naturally 
encompasses cross-saturation effects. 

For an unsaturated PMSM this energy function reads 

<t>q) = ^f^l + ^L~ q ^l 

where L d and L q are the motor self-inductances, and we 
recover the usual linear relations 



§d 



iq 9 2 Ul{(t)dAq) = J 1 ' 
L q 

Notice the expression for H should respect the symmetry 
of the PMSM w.r.t the direct axis, i.e. 

W(0d,-0g)=«(0d,0q), (5) 

which is obviously the case for Hi. Indeed ©-© is left 
unchanged by the transformation 

(Ud, U q , <j>d, <\> q , id, i q , u, 0, r L ) -> 

(u d , ~U q , <j) d , -<j) q , l d , ~iq, -UJ, ~9, ~T L ). 

C. Parametric description of magnetic saturation 

Magnetic saturation can be accounted for by considering 
a more complicated magnetic energy function H, having Hi 
for quadratic part but including also higher-order terms. From 
experiments saturation effects are well captured by considering 
only third- and fourth-order terms, hence 

H{(j)d,(t)q) = Hl((/)d,4>q) 

3 4 



i=0 



u q 1 

2=0 



This is a perturbative model where the higher-order terms 
appear as corrections of the dominant term Hi. The nine 
coefficients ol^ together with Ld, L q are motor dependent. 
But © implies a 2 ,\ = «o,3 = <^3,i = &i,3 — 0, so that the 
energy function eventually reads 

H((/>d, <f>q) = UMd, <t>q) + «3,O0rf + OLi i2 (t>d4> 2 q 

+ «4,o^ + a 2i2 (j) 2 d 4> 2 + «O,40g. (6) 
From © and © the currents are then explicitly given by 

id = P~ + 3a 3 ,o0d + a i^ 2 q + 4a 4 , o 0d + 2a 2 ,20d<^ (7) 



iq = Y~ + 2<*i,2<l)d<l>q + 2a 2 ,20d^9 + 4a 0,4^, 

which are the so-called flux-current magnetization curves. 



(8) 



To conclude, the model of the saturated PMSM is given 
by ©-© and ©-©, with <j) d ,<j) q ,uj,6 as state variables. 
The magnetic saturation effects are represented by the five 
parameters a^o, ^1,2, «4,o, ^2,2, «o,4- 

D. Model with id,i q as state variables 

The model of the PMSM is usually expressed with currents 
as state variables. This can be achieved here by time differen- 
tiating i dq =X dq {<t>dq), 



di d q _ nr u ,d<j>d q 



with given by ©. Fluxes are then expressed as <\> dq = 
^dqitdq) by inverting the nonlinear relations ©-©; rather 
than performing the exact inversion, we can take advantage 
of the fact the coefficients a^j are experimentally small. At 
first order w.r.t. the a^j we have <j)d = L d i d + @(\ a i,j\) anc * 
<j) q = L q i q + 0(\aij\); plugging these expressions into ©-© 
and neglecting (9(|o^j| 2 ) terms, we easily find 

<t>d = L d (id ~ 3a 3}0 L 2 d i 2 d - a 1}2 L 2 i 2 

- 4a 4j0 L 3 d i 3 d - 2a 2a L d L 2 q i d i 2 q ) (9) 

(j) q = L q (i q ~ 2(X\ , 2 L rfL q l d i q — 

2a 2 , 2 L 2 d L q i 2 d i q - 4a 0A L 3 q i 3 q ) . (10) 
Notice the matrix 

(%:<::! ?:<::>) <>■> 

with coefficients easily found to be 

Gdd(idq) = 7- + §a$$Ldid + ^2a^L 2 d i 2 d + 2a 2i2 L 2 i 2 

J^d 

G d q{idq) = 2ai y2 L q i q + ^Ct 2 ^ 2 L d i d^ qi q 

G qq (idq) = 7- + 2a\^ 2 L d i d + 2a 2 ,2^^d + 12a ,4^^, 
is by construction symmetric; indeed 

nT // \ _ / diiH(4>d, <j>q) d 21 H{(/)d,(t)q) 

dq[M ~\d 12 H^ d ,^q) d 22 H(^ d ,M 

and d\ 2 H = d 2 \H. Therefore the inductance matrix 

r L d d(idq) L dq {idq)\ . = (G d d{idq) G dq (id q )\ 1 

K Ld q (idq) L qq (i dq ) J \Gd q (idq) G qq (id q ) J 

is also symmetric, though this is not always acknowledged in 
saturation models encountered in the literature. 

III. Position estimation by high frequency voltage 

INJECTION 
A. Signal injection and averaging 

A general sensorless control law can be expressed as 

u a p = M e u l8 (12) 

(13) 



dOr. 



= a(MeJ lS ,0 c ,rj,t) 



dt 
dr\ 

~dt 



(14) 

(15) 
(16) 
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where the measured currents i a p = MqJ 7 § are used to 
compute u 7 §, uj c and the evolution of the internal (vector) 
variable r\ of the controller; C and uj c are known by design. 

It will be convenient to write the system equations ([T])-© 
in the 7 — S frame defined by x 7 s '= ^-Q-e c ^dq-> which gives 



D 1 5 



dt 

J duo 3 



u l5 - Ri l5 - oJcJC^s - uJCMe-e^m (17) 



n 2 dt ~ i^ 5 ^ 15 + M °-°c^ 
dO 

Tt=^ 



(18) 
(19) 



where from ©-([8]) currents and fluxes are related by 

V = M e _ 9( X dq (Mj_ e ^ l5 ). (20) 

To estimate the position we will superimpose on some 
desirable control law (IT6b a fast-varying pulsating voltage, 



(21) 



where / is a 27r-periodic function with zero mean and u^s 
could like U 7 s depend on Mo c i 1 5,6 c ,ri,t (though it is al- 
ways taken constant in the sequel). The constant pulsation 
Q is chosen "large", so that f(ftt) can be seen as a "fast" 
oscillation; typically ft := 2tt x 500rad/s in the experiments 
in section HV-Cl 

If we apply this modified control law to (fT71)-([T9l), then it 
can be shown the solution of the closed loop system is 







(22) 


CO = UJ -f 




(23) 


(9 = 5 + 




(24) 


#c = 6c H 




(25) 


?7 = 77 + 




(26) 



where F is the primitive of / with zero mean (F clearly has 
the same period as /); (<j> s ,uJ, 0, C , rj) is the "slowly- varying" 
component of (0 7< y, 0, 6 c ,rj), i.e. satisfies 



d(\) 



7(5 



dt 



= u 7 s — Ri^s — u c JC(j)^ § — ujJCM-q_q 



J duo 3-t w — 
~dt_ 
dO 



where 



—-=uj 
dt 

dO c 

dn , - - _ , 

— = a{M^ c i 76 ,O c ,r] J t), 



lu c = Q JC (M-Q c i l5 ,Q c ,rj,t) 
u l5 =Uys(M0 i 7 s,O c ,V>t)' 



(27) 



Notice this slowly-varying system is exactly the same as (TTTb— 
(|T9t acted upon by the unmodified control law (fT2l)-([T6l). In 
other words adding signal injection: 

• has a very small effect of order 0{-^) on the mechanical 
variables Q,uu and the controller variables Cl r] 

• has a small effect of order O(^) on the flux <p 1 5\ this 
effect will be used in the next section to extract the 
position information from the measured currents. 

The proof relies on a direct application of second-order 
averaging of differential equations, see l37ll section 2.9.1 
and for the slow-time dependance section 3.3. Indeed setting 

£ := h> a := f' anc * x := OAytf* ^' v)> (ITTb— (|T9b acted 
upon by the modified control law (IT2l)-(IT31) and (f2~Tl) is in 
the so-called standard form for averaging (with slow-time 
dependance) 

— = efi(x,ea,a) := e(f 1 (x,ea)-\-f 1 (x,ea)f(a)), 

with fi T-periodic w.r.t. its third variable (T = 2tt in our 
case) and e as a small parameter. Therefore its solution can 
be approximated as 

x(cr) = z(a) + s(ui(z(a),scr, a) + 0(e 2 ), 

where z(a) is the solution of 

dz 



and 



1 f T 

gi (y, ecr) := — J f x (y, eg, s)ds = f 1 (y, £<r) 

vxiy.ecr.cr) := / (f 1 (y,sa,s)-g 1 (y,scr))ds 
Jo 

= fi(y,£<r) / /(s)ds 



1 r T 

u^y.ea.a) := v 1 (y : ea,a) - — / v 1 (y, ecr, s)ds 

1 Jo 

= fi(y,ea)F(a) 
K 2 (y, server) := difi(y,ea,a)ui(y,ea,a) 
- diU!(y,ea, o)gi{y,eo) 
= UiJi](y,ea)F(a) 

+ ^i/i(y,£(j)/i(i/,g(7) — 

1 r T 

^2(2/, ecr) := - y ^2(2/, ecr, 5)^8 = 0. 
We have set 

[7l 5 /l] (y , ^) : = 3l7l (2/ j (2/ j ^) - ^1 h (2/, ^)7! (2/, SCr) 

and F(a) := f(s)ds — ^ J Q T f(s)dsda, i.e. F is the (of 
course T-periodic) primitive of / with zero mean. 

Translating back to the original variables eventually yields 
the desired result (I22M261). 
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B. Position estimation 

We now express the effect of signal injection on the cur- 
rents: plugging (l22l) into (l20l) we have 



i l5 - M e _ 0c+o ^ ) 

X dq[ 



1 



V + * 74 F(fit) + 0(^), 



(28) 



where we have used (1771) and performed a first-order expansion 
to get 

I 7 , := M- e _- e m dq (Mf_^ s )Ml- dc ^- 



Ms 



T "v)K-*.if- (29) 



We will see in the next section how to recover i 7 s and i 7( $ 
from the measured currents z 7 £. Therefore (l29l) gives two 
(redundant) relations relating the unknown angle to the 
known variables C , z^, z 7 £, S^g, provided the matrix 

i 7 i) := M^DIrf, (I" 1 (Mji 7 ,)) Mj 

effectively depends on its first argument \i. This "saliency con- 
dition" is what is needed to ensure nonlinear observability. The 
explicit expression for <S(/i, i 7 §) is obtained thanks to ([TH . In 
the case of an unsaturated magnetic circuit this matrix boils 
down to 




Mi 



^^cos2/i 

L d ~\-L q r 



and does not depend on \ notice this matrix does not depend 
on /i for an unsaturated machine with no geometric saliency. 
Notice also (l29t defines in that case two solutions on ] — 7r, it] 
for the angle 6 since S(fi,i y s) is actually a function of 2/i; 
in the saturated case there is generically only one solution, 
except for some particular values of i 7 s- 

There are several ways to extract the rotor angle information 
from (|29K especially for real-time use inside a feedback law. In 
this paper we just want to demonstrate the validity of d29l) and 
we will be content with directly solving it through a nonlinear 
least square problem; in other words we estimate the rotor 
position as 



arg mm 

MG]— 7T,7 



C. Current demodulation 



i 1 5 - S(fj,,i 7 s) 



Uj5_ 

7T 



(30) 



To estimate the position information using e.g. (l30l) it is 
necessary^ to extract the low- and high-frequency components 
i 1 s and i 1 $ from the measured current i^. Since by (f28l) 

z 7 s(t) « i 1 5(t) J ri 1 s(t)F(Vtt) with i 7 £ and i 7 ^ by construction 



nearly constant on one period of F, we may write 

1 i l 

~ 7^ / i<ys(s)ds 

1 Jt-T 



fZF*(Sls)ds 



where T := ^y. Indeed as F is 27r-periodic with zero mean, 

i 1 s{s)ds ^ i 1 s(t) I ds J ci 1 $(t) I F(tts)ds 

t-T Jt-T Jt-T 

= Ti lS (t) 



I i 7 s(s)F(ns)ds&i 7 s(t) [ F(Qs)ds 

Jt-T Jt-T 

+ ! 7 s(t) [ F 2 (ns)ds 

Jt-T 

= l yS (t) [ F 2 (ns)ds. 
Jo 
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Fig. 1. Experimental time response of id in f3H-f32l 

IV. Estimation of magnetic parameters 

The seven parameters in the saturation model ©-([8]) must 
of course be estimated. This can be done with a rather simple 
procedure also relying on signal injection and averaging. 

A. Principle 

The rotor is locked in the position 6 := 0, hence the model 
(Q])-© reduces to uj = and 



dt 



= u dq — Ridq, 



(31) 



with id q = Tdq{(j>dq)' Moreover Ud q can now be physically 
impressed and id q physically measured. 

As in section IIII-A1 but now working directly in the d — q 
frame, we inject a fast- varying pulsating voltage 

Udq = U dq + Udq f (tot), (32) 

with constant Ttd q and Udq- The solution of d3Tl)-(l32b is then 



TABLE I 

Rated and estimated magnetic parameters of test motors 



with u q = 0, hence i q = 0, and u q — (l34t reads 



Motor 


IPM 


SPM 


Rated power 


750 W 


1500 W 


Rated current I n (peak) 


4.51 A 


5.19 A 


Rated voltage (peak per phase) 


110 V 


245 V 


Rated speed 


1800 rpm 


3000 rpm 


Rated torque 


3.98 Nm 


6.06 Nm 


n 


3 


5 


R 


1.52 Q 


2.1 O 


A (peak) 


196 mWb 


155 mWb 


L d 


9.15 mH 


7.86 mH 


L q 


13.58 mH 


8.18 mH 




0.039 


0.056 


OLl^LdLqln 


0.053 


0.055 




0.0051 


0.0164 




0.0171 


0.027 


a 0A Lpl 


0.0060 


0.0067 



where (j) d , the "slowly- varying" component of (j>d q , satisfies 



dt 



u d - Rid, 



with id q = Xd q ((j)dq)- Moreover ([29b now boils down to 



i dq = Dl dq (l dq \id q ))^ L . 



n 



(33) 



(34) 



Since Ud q is constant (l33t implies Ridq tends to Ttd q , hence 
after an initial transient id q is constant. As a consequence id q 
is by (l34h also constant. Fig. [T] shows for instance the time 
response of id for the SPM motor of section IIV-CI starting 
from id(0) = and using a square function /; notice the 
current ripples seen on the scope are max rG [o 5 27r] ^X r ) = f 
(since / is square with period 2tt) smaller than id q . 

The magnetic parameters can then be estimated repeatedly 
using (f34l) with various values of Ud q and Ud q , as detailed in 
the next section. 

B. Estimation of the parameters 

From (ITTb the entries of DXd q (T^(id q )) are given by 



Gdd(idq) 



1 

L~ d 



-2 



• 6a3,oL d i d + l2a±$L d i d + 2a 2 , 2 L q i q 



Gd q (idq) = 2a 1)2 L q iq + 4a 2y2 L d idL q i q 
- 1 

G qq (idq) = — + 2a lj2 L d i d + 2a 2 ,2^d*d + 12 «0,4^^q 

Since combinations of the magnetic parameters always enter 
linearly those equations, they can be estimated by simple linear 
least squares; moreover by suitably choosing Ud q and Ud q , the 
whole least squares problem for the seven parameters can be 
split into several subproblems involving fewer parameters: 
• with Udq := 0, hence idq — 0, (l34t reads 

1 u d 



L d = 



id 



T _ 1 U Q 



(35) 
(36) 



V± ( — 



i q =0 



with = 0, hence id := 0, and u q = (l34l) reads 

2d / 1 



^d 



2u d 




-150 -100 -50 50 100 150 200 

% of rated current id 




-200 



-100 100 

% of rated current i q 



200 



Fig. 2. IPM: fitted values vs measurements for d37l and d39l 
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Fig. 3. SPM: fitted values vs measurements for J37) and (39) 
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with Ud 



0, hence id 

2lt n 



0, and Ud = (l34t reads 



*d -FTOLxoLqiq 



12(^0,4^^0 



(40) 
(41) 



and L q are then immediately determined from (f35l) 
and (f36b ; o^o and a^o are jointly estimated by least squares 
from d37l) ; 0:2,2* ^1,2 and ao,4 are separately estimated by least 
squares from respectively ([38]), d39])-(|40|) and (I4T1) . 

C. Experimental setup 

The methodology developed in the paper was tested on 
two types of motors, an Interior Magnet PMSM (IPM) and a 
Surface-Mounted PMSM (SPM), with rated parameters listed 
in the top part of table U 

The experimental setup consists of an industrial inverter 
(400 V DC bus, 4 kHz PWM frequency), an incremental enco- 
der, a dSpace fast prototyping system with 3 boards (DS1005, 
DS5202 and EV1048), and a host PC. The measurements 



are sampled also at 4 kHz, and synchronized with the PWM 
frequency. The load torque is created by a 4 kW DC motor. 



D. Estimation of the magnetic parameters 

We follow the procedure described in section [IVJ with the 
rotor locked in the position := 0, a square wave voltage 
with frequency Vt := 2tt x 500rad/s and constant amplitude 
u d or u q (15 V for the IPM, 14 V for the SPM) is applied 
to the motor; but for the determination of Ld,L q where 
Ud = u q := 0, several runs are performed with various 
Ud (resp. u q ) such that id (resp. i q ) ranges from —200% to 
+200% of the rated current. The magnetic parameters are then 
estimated by linear least squares according to section IIV-BI 
yielding the values in the bottom part of table H Notice the 
SPM exhibits as expected little geometric saliency (Ld ~ L q ) 
hence the saturation-induced saliency is paramount to estimate 
the rotor position. Notice also the cross-saturation term ol\i is 
as expected quantitatively important for both motors. 

The good agreement between the fitted curves and the 
measurements is demonstrated for instance for (l37l) and 










r \ 













20 40 60 80 100 120 140 160 180 200 



(c) 



^ 100 



20 40 60 80 100 120 140 160 180 200 




20 40 60 80 100 120 140 160 180 200 
time in sec 



20 40 60 80 100 120 140 160 180 200 
time in sec 



Fig. 4. Long test under various conditions for IPM: (a) measured — C , 
estimated 6 — 9 C with and without saturation model; (b) measured speed uo, 
reference speed oo c ; (c) load torque r^; (d) voltages u 7 ^ 



Fig. 5. Long test under various conditions for SPM: (a) measured — 9 C , 
estimated — C with and without saturation model; (b) measured speed uo, 
reference speed uj c ; (c) load torque r^; (d) voltages u^i 
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on Fig. [JEl notice (f3Tb corresponds to saturation on a single 
axis while (f39l) corresponds to cross-saturation. 

E. Validation of the rotor position estimation procedure 

The relevance of the position estimation methodology de- 
veloped in section Jill is now illustrated on the two test motors, 
using the parameters estimated in the previous section. Since 
the goal is only to test the validity of the angle estimation 
procedure, a very simple V/f open-loop (i.e. Q c and U 7 s do 
not depend on i 7 §) control law is used for (IT2T)— (IT6b ; a fast- 
varying (ft := 2tt x 500rad/s) square voltage with constant 
amplitude is added in accordance with (Ell) , resulting in 

-dt =Uc{t) 

u l5 = u r * s (t) + cj c (t)0 m + u l5 f(Vtt). 

Here oj c (t) is the motor speed reference; u^(t) is a filtered 
piece-wise constant vector compensating the resistive voltage 
drop in order to maintain the torque level and the motor 
stability; finally u 7 $ := (2, 0) T with u := 15 V. 

The rotor position is then estimated according to (l30l) . 



1 ) Long test under various conditions, Fig. @[3} Speed and 
torque are changed over a period of 210 seconds; the speed 
remains between ±5% of the rated speed and the torque varies 
from 0% to 180% of the rated toque. This represents typical 
operation conditions at low speed. 

When the saturation model is used the agreement between 
the estimated position and the measured position is very 
good, with an error always smaller than a few (electrical) 
degrees. By contrast the estimated error without using the sat- 
uration model (i.e. with all the magnetic saturation parameters 
aij taken to zero) can reach up to 40° for the IPM and 70° 
the SPM. This demonstrates the importance of considering 
an adequate saturation model including in particular cross- 
saturation. 

2 ) Slow speed reversal, Fig. \fi§7j[ This is an excerpt of the 
long experiment between 35 s and 55 s. The speed is slowly 
changed from -0.2% to +0.2% of the rated speed at 150% of 
the rated torque. This is a very demanding test since the motor 
always remains in the poor observability region, moreover 
under high load. Once again the estimated angle closely agrees 
with the measured angle. 
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Fig. 6. Slow speed reversal for IPM: (a) measured 0, estimated 0; (b) Fig. 7. Slow speed reversal for SPM: (a) measured 6, estimated 6; (b) 
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3 ) Load step at zero speed, Fig. [SHPf This is an excerpt of 
the long experiment around t = 125 s. The load is suddenly 
changed from 0% to 100% of the rated torque while the motor 
is at rest. This test illustrates the quality of the estimation also 
under dynamic conditions. 

V. Conclusion 

We have presented a simple parametric model of the satu- 
rated PMSM together with a new procedure based on signal 
injection for estimating the rotor angle at low speed relying on 
an original analysis based on second-order averaging. This is 
not an easy problem in view of the observability degeneracy 
at zero speed. The method is general in the sense it can 
accommodate virtually any control law, saturation model, and 
form of injected signal. The relevance of the method and the 
importance of using an adequate magnetic saturation model 
has been experimentally demonstrated on a SPM motor with 
little geometric saliency as well as on an IPM motor. 
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